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Abstract 

Experimental and theoretical progress concerning the rare decay 7r° — > e + e~ is briefly reviewed. 
It includes the latest data from KTeV and a new model independent estimate of the decay branching 
which show the deviation between experiment and theory at the level of 3.3o\ The predictions for 
n and 7/ decays into lepton pair are presented. We also comment on the impact on the pion rare 
decay estimate of the data of BABAR collaboration on the pion transition form factor at large 
momentum transfer. 

PACS numbers: 
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I. INTRODUCTION 



Astrophysics observables tell us that 95% of the matter in the Universe is not described 
in terms of the Standard Model (SM) matter. Thus, the search for the traces of New Physics 
is a fundamental problem of particle physics. There are two strategies to look for the effects 
of New Physics: experiments at high energy and experiments at low energy. In high-energy 
experiments it is considered that due to a huge amount of energy the heavy degrees of 
freedom presumably characteristic of the SM extension sector are possible to excite. In low- 
energy experiments it is huge statistics that compensates the lack of energy by measuring 
the rare processes characteristic of such extensions. At present, there is no any evidence for 
deviation of SM predictions from the results of high-energy experiments and we are waiting 
for the LHC epoch. On the other hand, in low-energy experiments there are rough edges 
indicating such deviations. The most famous example is the muon (g — 2). Below it will be 
shown that due to recent experimental and theoretical progress the rare process tt° — > e + e~ 
became a good SM test process and that at the moment there is a discrepancy between the 
SM prediction and experiment at the level of 3.3<r deviation. 



II. KTEV DATA 

In 2007, the KTeV collaboration published the result lj for the branching ratio of the 
pion decay into an electron-positron pair 

B™L (vr - e+e-) = (7.48 ± 0.38) • lO" 8 . (1) 

The result is based on observation of 794 candidate 7r° — > e + e~ events using Kl — > 3tt° as 
a source of tagged 7r°s. Due to a complicated chain of the process and a good technique for 
final state resolution used by KTeV this is a process with low background. 



III. CLASSICAL THEORY OF tt° -> e + e DECAY 

The rare decay T » _ has b tudied theoretically over the years , starting Wlth 

the first prediction of the rate by Drell |2J. Since no spinless current coupling of quarks to 
leptons exists, the decay is described in the lowest order of QED as a one-loop process via the 
two-photon intermediate state, as shown in Fig. 1. A factor of 2 (m e /m n ) 2 corresponding 
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to the approximate helicity conservation of the interaction and two orders of a suppress the 
decay with respect to the tt° — > 77 decay, leading to an expected branching ratio of about 
10~ 7 . In the Standard Model contributions from the weak interaction to this process are 
many orders of magnitude smaller and can be neglected. 



Ai) 




FIG. 1: Triangle diagram for the 7r° — > e+e process with a pion tt° — > 7*7* form factor in the 
vertex. 



To the lowest order in QED the normalized branching ratio is given by 



R _> e+e") 



B (n° -> e+e-) ' x :> 



2(-^) P e (ml)\A(ml) 
B (tt° ^77) V 7T m, ' v *' 1 K J 



(2) 



where /3 e (q 2 



1 - 4%, 5 (tt° -> 77) = 0.988. The amplitude A can be written as 



{ 2\ f ^ ^ 7~i 



(fc 2 ,(A;-g) 2 ) 



g 2 fc 2 - (gfc) 5 



(A; 2 + is) ((k — q) + is) ((k — p) — m 2 e +is) 



, (3) 
(4) 



where q 2 = m 2 ,p 2 = m 2 . -F 7r7 * 7 * is the form factor of the transition tt° — > 7*7* with off-shell 
photons. 

The imaginary part of A is defined uniquely as 



Im.A (q 2 ) 
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2/3 e (q 2 



In {y e (g 2 )) 



(5) 



1 - (3 e (q 2 ) 

ye{q) = TTKWy 



It comes from the contribution of real photons in the intermediate state and is model inde- 
pendent since -F 7r7 » 7 » (0, 0) = 1. Using inequality \A\ 2 > (Inu4.) 2 one can get the well-known 
unitary bound for the branching ratio 



B (tt° -> e+e') > B unitmy (tt° -> e + e~) = 4.69 • 1(T 8 . 



(6) 



One can attempt to reconstruct the full amplitude by using a once-subtracted dispersion 
relation [s| 

If one assumes that Eq. (J5j) is valid for any q 2 , then one arrives for q 2 > Am 2 at p, |7|, |8[] 



Re^ (q 2 ) =A(q 2 = 0) + 



Pe(q 2 



\ln 2 (y e (q 2 )) + ^ + U 2 (-y e (q 2 )) 



(8) 



where Li 2 (z) = — J* (dt/t) In (1 — t) is the dilogarithm function. The second term in Eq. (JHJ) 
takes into account a strong q 2 dependence of the amplitude around the point q 2 = occurring 
due to the branch cut coming from the two-photon intermediate state. In the leading order 
in (m e /m n ) 2 , Eq. (JSj) reduces to 



HMK)^ (*» = <>)+ 



(9) 



Thus, the amplitude is fully reconstructed up to a subtraction constant. Usually, this 
constant containing the nontrivial dynamics of the process is calculated within different 
models describing the form factor F^k 2 ^ 2 ) 



shown in 



lOj. However, it has recently been 



10( that this constant may be expressed in terms of the inverse moment of the 



pion transition form factor given in symmetric kinematics of spacelike photons 



A (q 2 = 0) = 3 In 



5 _ 3 
4 ~ 2 



M TP 



+ 



dt 



t 



(10) 



Here, /i is an arbitrary (factorization) scale. One has to note that the logarithmic dependence 
of the first term on [i is compensated by the scale dependence of the integrals in the brackets. 
In this way two independent processes becomes related. 



IV. IMPORTANCE OF CLEO DATA ON F 7r7 * 7 



In order to estimate the integral in Eq. (flUj) . one needs to define the pion transition form 
factor in symmetric kinematics for spacelike photon momenta. Since it is unknown from 
the first principles, we will adapt the available experimental data to perform such estimates. 
Let us first use the fact that F,r 7 * 7 * (t, t) < -F 7r7 . 7 . (t, 0) for t > in order to obtain the 
lower bound of the integral in Eq. f fTOl) . For this purpose, we take the experimental results 
from the CELLO ll| and CLEO 12| Collaborations for the pion transition form factor in 
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asymmetric kinematics for spacelike photon momentum which is well parameterized by the 
monopole form [3] 



-^7*7* (t, 0) — ,/ CLEO' (11) 



1 

s CLEO = ( 77 g ± 22 MeV )2 



For this type of the form factor one finds from Eq. (fTUj) that 

o / „CLEO \ c 

-4(? 2 = 0) >-|ln(^^a- r j-| = -23.2 ± 0.1. (12) 

Thus, for the branching ratio we are able to establish the important lower bound which 
considerably improves the unitary bound given by Eq. (J6j) 

B (vr _> e + e") > 5 CLEO (tt° -> e+e") = (5.84 ± 0.02) ■ lO" 8 . (13) 

It is natural to assume that the monopole form is also a good parametrization for the 
form factor in symmetric kinematics 

F, 7 . 7 .(M) = TT ^ 7 . (M) 

The scale s\ can be fixed from the relation for the slopes of the form factors in symmetric 
and asymmetric kinematics at low t [3], 

dF^try* (t, t) I 2 ^^ 7T7*7* 0) 



dt \ t=0 - dt 



(15) 

t=o 

that gives Si = so/2. Note that a similar reduction of the scale is also predicted by OPE 
QCD from the large momentum behavior of the form factors: Si PE = Sq PE /3 [l4j. Thus, 
the estimate for A (0) can be obtained from Eq. (fl2j) by shifting the lower bound by a 
positive number which belongs to the interval [3 ln(2)/2, 3 ln(3)/2] 



(16) 



With this result the branching ratio becomes 

B (tt° -> e+e") = (6.23 ± 0.09) • 10~ 8 . (17) 
This is 3.3 standard deviations lower than the KTeV result given by Eq. (JT]). 
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V. OTHER DECAY MODES 



The r] — > / + Z decay can be analyzed in a similar manner. As in the pion case, the CLEO 



Collaboration has parameterized the data for the rj- meson in the monopole form 12] : 



cleo 1 
Frr/*7* ^) = ]_ _|_ £/ s CLEO' 

S CLEO = ( 7?4 ± 2g Me y^2 ^ 

which is very close to the relevant pion parameter. Then following the previous case (with 
evident substitutions), one finds the bounds for the q 2 — > limit of the amplitude rj — > fi + fi~ 

as 

o / s CLEO \ r 

A if = 0) > ~\ In M^-J - | = - (7.2 ± 0.1) , (19) 

and for 77 — > e + e~ one gets again Eq. ( fl2i) . The obtained estimates allow one to find the 
bounds for the branching ratios 

B (rj -> < (6.23 ±0.12) ■ 10" 6 , (20) 

5 (rj -> e + e~) > (4.33 ± 0.02) • 10" 9 . 

It is important to note that for the decay rj — > A i+ /^~ we get the upper limit for the branching. 
This is because the real part of the amplitude for this process taken at the physical point 
q 2 = m 2 for the parameter Sq^ eo remains negative and a positive shift due to the change of 
the scale so v — > si v reduces the absolute value of the real part of the amplitude |Re^4. [mtj |. 
At the same time, considering the decays of 7r° and rj into an electron-positron pair, the 
evolution to physical point (jHD makes the real part of the amplitude to be positive for 
the parameter So LEO and the absolute value of the real part of the amplitude increases in 
changing the scales of the meson form factors. 

The predicted branchings are given in the Table and compared with existing experimental 
data. The so-called unitary bound appears if in (j2J) only the imaginary part of the amplitude 
which is model independent is taken into account. The CLEO bound corresponds to the 
estimate of the real part of the amplitude basing on the CELLO and CLEO data on the 
meson transition from factors [10]. The fourth column of the Table contains the predictions 



where in addition the constraint from O 
photon virtualities is taken into account 



3 E QCD on the transition form factor for arbitrary 
3- 



The results that take into account the mass corrections to the amplitude [20] , mainly due 
to powers (M/A) 2 , where M is the pseudoscalar meson mass and A « M p is characteristic 
scale of the pion transition form factor, are given in the fifth column of the Table. They 
are essentially visible for rj (77') meson decays. It is interesting that for rj decay to muons 
the mass correction shifts the theoretical prediction in the direction to the unitary bound 
and thus opposite to the experimental result [lfj]. Thus, it would be very desirable to 
check experimentally the predicted bounds for the process rj —>■ yU + /i~. Also note the recent 
measurement by the WASA/Celcius Collaboration [jjj] which improves the upper limit for 
the branching rj — > e + e~. 

For rf decays there appear new thresholds in addition to the two-photon one. In general, 
this violates the unitary bound because the correction A z gains an imaginary part at z = 
(M/A) 2 > 1 

+ 



In 



e(z-i). 



(21) 



As it is seen from the table, this happens for the 7/ — > fi/i channel. Nevertheless, it turns 
out that the predictions for this channel are quite accurate. 



TABLE I: Values of the branchings B (P — > l + l ) obtained in our approach and compared with 
the available experimental results. 



Rq 


Unitary bound 


CLEO bound 


CLEO±OPE [10] 


W 


Experiment 


R (vr° -► e+e") x 10 8 


> 4.69 


> 5.85 ±0.03 


6.23 ±0.12 


6.26 


7.49 ±0.38 [1] 


Ro (v -» M"V~) x 10 6 


> 4.36 


< 6.23 ±0.12 


5.12 ±0.27 


4.64 


5.8 ±0.8 [16] 


R (77 -» e+e-) x 10 9 


> 1.78 


> 4.33 ± 0.02 


4.60 ± 0.09 


5.24 


< 2.7- 10 4 [15J 


R (if -> ^ + /i~) x 10 7 


> 1.35 


< 1.44 ±0.01 


1.364 ±0.010 


1.30 




^(77' -» e+e~) x 10 10 


> 0.36 


> 1.121 ±0.004 


1.178 ±0.014 


1.86 





VI. POSSIBLE EXPLANATIONS OF THE EFFECT 

Therefore, it is extremely important to trace possible sources of the discrepancy between 
the KTeV experiment and theory. There are a few possibilities: (1) problems with (statistic) 
experiment procession, (2) inclusion of QED radiation corrections by KTeV is wrong, (3) 
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4) effects of new physics. At the mo- 



ll, the contribution of QED radiative 



unaccounted mass corrections are important, and 
ment, the last possibilities were reinvestigated. In 
corrections to the ir° — > e + e~ decay, which must be taken into account when comparing the 
theoretical prediction (1171) with the experimental result ([I]), was revised. Comparing with 
earlier calculations [18|, the main progress is in the detailed consideration of the 7*7* — > e + e~ 
subprocess and revealing of dynamics of large and small distances. Occasionally, this num- 
ber agrees well with the earlier prediction based on calculations [3] and, thus, the KTeV 



analysis of radiative corrections is confirmed. In |19l . |20] it was shown that the mass correc- 
tions are under control and do not resolve the problem. So our main conclusion is that the 
inclusion of radiative and mass corrections is unable to reduce the discrepancy between the 
theoretical prediction for the decay rate (fP7j) and experimental result ([T]). 

VII. 7T° -» e + e~ DECAY AS A FILTERING PROCESS FOR LOW MASS DARK 
MATTER 



If one thinks about an extension of the Standard Model in terms of heavy, of an order of 
100 GeV or higher, particles, then the contribution of this sort of particles to the pion decay 
is negligible. However, there is a class of models for description of Dark Matter with a mass 
spectrum of particles of an order of 10 MeV 26| . This model postulates a neutral scalar dark 
matter particle \ which annihilates to produce electron/positron pairs: XX ~* e + e~. The 
excess positrons produced in this annihilation reaction could be responsible for the bright 
511 keV line emanating from the center of the galaxy 27j. The effects of low mass vector 



28j | where the excess of 



boson U appearing in such model of dark matter were considered in 
KTeV data over theory put the constraint on coupling which is consistent with that coming 
from the muon anomalous magnetic moment and relic radiation [29]. Thus, the pion decay 
might be a filtering process for light dark matter particles. 



VIII. CONSTRAINTS FROM BABAR DATA 



Let us now consider the amplitude 7r° — > e + e in the context of the constituent quark 



model 



22 



23] . Within this model, the pion form factor is given by the quark-loop (triangle) 



(Fig. 1) diagram with momentum independent constituent quark mass. The result for the 
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form factor in symmetric kinematics is given by 23] 

Pat 



2M2, 3 n + 1 

F^^^-^ln^l-, (22) 



where (3n = \ll + ^y 2 -, and the amplitude at zero is 



It is interesting to remind that in asymmetric kinematics the form factor has double loga- 



rithmic asymptotics at large momentum transfer log 2 (Q 2 /M 2 ) 23J 



Recently BABAR collaboration announced new data [2l|] on pion transition form factor 
F V11 * (t, 0) at very large t up to 40 GeV 2 . The data are in strong contradiction with massless 
QCD predicting asymptotic 1/t behavior with single logariphmic modulation in leading 
order. However, the BABAR data are well fitted by ( |24|) if the parameter Mq is taken 
as Mq = 135 MeV. Then the amplitude at zero is .4.(0) = —20.98 which is still far from 
numbers corresponding to KTeV: A (0) = —18.6 ±0.9. In [24J it was also noted that in order 
to get reasonable estimate of the hadronic vacuum polarization to the anomalous magnetic 
moment of muon it is necessary to take the quark mass rather small Mq = 180 MeV, the 
value which is close to the pion mass. 

Within perturbative QCD one possible scenario to explain the form factor growth faster 
than expected is to assume that the full resummation of perturbative logarithms is impor- 



tant 1 . Indeed, sometime ago Manohar augmented 25j that for the pion transition form 
factor in kinematics when one photon is virtual and other is quasireal the resummation ef- 
fects dominate its behavior at large fixed Q 2 and probably BABAR (if confirmed) observes 
these effects (logarithmic or powerlike). 
5 

Further independent experiments at KLOE, NA48, WASAatCOSY, BES III and other 
facilities will be crucial for resolution of the problem with the rare leptonic decays of light 
pseudoscalr mesons. Also it is important to confirm the theoretical base for maximally model 



1 I am indebted to A. A. Pivovarov for discussion this point. 
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n 

FIG. 2: The transition form factor vr° -> 7*7. The data are from CELLO 11], CLEO [12] and 
BABAR |2l| Collaborations. The dashed line is massless QCD asymptotic limit. The solid line is 
the form factor calculated from massive triangle diagram, Eq. (|24p . with mass parameter taken as 
M q = 135 MeV. 

independent prediction of the branchings by getting more precise data on the pion transition 
form factor in asymmetric as well in symmetric kinematics in wider region of momentum 
transfer that is soon expected from BABAR and BELLE collaborations. 
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